Gauge invariance of the dual Meissner effect in QCD 
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The dual Meissner effect is described and numerically observed in a gauge- invariant way in lattice 
Monte-Carlo simulations of pure SU(2) QCD. A gauge-invariant Abelian-like field strength is defined 
in terms of a unit- vector in color space which is constructed by a non-Abelian field strength itself. A 
gauge-invariant monopole-like quantity is defined by a violation of the Bianchi identity with respect 
to the Abelian-like field strength. The squeezing of the non-Abelian electric field \/^2 a (Ef) 2 between 
a pair of static quark and anti-quark occurs due to the solenoidal current coming from the gauge- 
invariant monopole-like quantity. An equation similar to the dual London equation is confirmed 
approximately in the long-range region. 
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One of the most essential problems of color confine- 
ment in QCD is to explain the mechanism of the flux 
squeezing of non-Abelian electric fields between a pair of 
static quark and anti-quark. In SU(2) QCD, ^2 a (Ef) 2 or 
y/^2 a (Ef ) 2 is expected to be squeezed to reproduce the 
linear static potential. Numerically the expected squeez- 
ing of the gauge-invariant combination of the electric field 
was observed beautifully in lattice ST/ (2) QCD 0. 

Thirty years ago, 'tHooft and Mandelstam js] con- 
jectured that the dual Meissner effect is the color con- 
finement mechanism of QCD. However what causes the 
dual Meissner effect and how to treat the non-Abelian 
property were not clarified. 

An interesting idea is to utilize a topological monopole 
like the 'tHooft-Polyakov monopole 0, || . An important 
quantity is a 'tHooft field strength f^ v — J2 a + 
J2abc e abc'n a (D fl n) b (D v n) c where n a is a unit vector com- 
posed of gluonic fields transforming as an adjoint repre- 
sentation in color space and (D^n) b is a covariant deriva- 
tive. A monopole picture can be seen more clearly if we 
project SC/(3) QCD to an Abelian U{1) 2 theory by a par- 
tial gauge fixing 6j. Then we have an Abelian U(l) 2 the- 
ory with Abelian electric and magnetic charges. It is con- 
jectured in Ref.0] that the condensation of the Abelian 
monopoles causes the dual Meissner effect explaining the 
color confinement. 

However there is a serious problem in this scenario. 
Namely there exist infinite ways of choosing n a or in other 
words infinite possible Abelian projections. Moreover, 
the monopole condensation, if happens, can explain only 
the squeezing of an Abelian-like electric field fa. How 
good an approximation it is to the real and expected flux 
squeezing of \/^Z a {Ef) 2 depends strongly on the choice 
of 

An Abelian projection adopting a special gauge called 
Maximally Abelian gauge (MA) 0, HEj is found to give 
us interesting results [lCt llll Il2| supporting importance 
of the Abelian monopoles. In this case, the Abelian elec- 
tric field can approximate very well the long-range behav- 
ior of the non-Abelian one, since off-diagonal components 



are suppressed. However such beautiful results are not 
seen in other general gauges. 

It is the purpose of this note to show numerically 
that the dual Meissner effect is observed in a gauge- 
invariant way with the use of a gauge-invariant Abelian- 
like field strength and a monopole-like quantity. We 
do not need any Abelian projection nor any gauge- 
fixing. Monte-Carlo simulations of quenched SU (2) QCD 
are performed. It is found that the squeezing of the 
non-Abelian electric field \f}2 a {Ef) 2 occurs and the 
solenoidal current from the gauge-invariant monopole- 
like quantity is responsible for the flux squeezing. The 
magnetic displacement current observed previously in 
Landau gauge |l3j is found to be negligible. Preliminary 
results are obtained with respect to the vacuum type of 
the confinement phase. The SU(2) QCD vacuum seems 
near the border between the type 1 and the type 2 dual 
superconductors. The present numerical results are not 
perfect, since the continuum limit, the infinite- volume 
limit and the real St/(3) case are not studied yet. Nev- 
ertheless the authors think the results obtained here are 
very interesting to general readers, since they show for 
the first time the flux squeezing of non-Abelian electric 
fields is working in a gauge-invariant way due to the dual 
Meissner effect without performing any Abelian projec- 
tion. 

Let us define an Abelian-like field strength and a 
gauge-invariant monopole-like quantity in QCD. The 
field strength is written in terms of a unit-vector in 
color space which is constructed by a non-Abelian field 
strength itself: 



\n a (x)F a (x). 



(1) 



n." is a unit vector 



F^ v is a non-Abelian field strength, 
in color space transforming as an adjoint representation 
m su{2) OH E3 and is a symmetric tensor in space 



i° v (x) = e^F^x)/ r£(F>„(x)) 2 , 



(2) 



where e^ v is an antisymmetric tensor with a sign con- 
vention e fl<l/ = 1. The opposite sign convention can be 
adopted. Note that Eq.© is just equal to the gauge- 
invariant absolute value of the non-Abelian field strength 
itself except for the sign. Hence it is not a simple Lorenz 
tensor. It is noted that an electric field component Ei 
defined by fa is -^EaK) 2 tne squeezing of which is 
to be explained. 

A gauge-invariant monopole-like quantity is defined by 



1 



(3) 



which is conserved but is not a simple Lorenz vector. 
Hereafter we call the monopole-like quantity simply as 
'monopole'. We get from Eq.© 



VxEl d A B = Airk, 

V B = -47rfc 4 , 



(4) 
(5) 



where Bi = 1/2 ey/./^ and the vector notation is 
with respect to the three-dimensional space. Note that 
the magnetic charge defined in Eq. © does not satisfy the 
Dirac quantization condition with respect to bare charges 
contrary to the usual case of a magnetic charge defined 
in terms of a 'tHooft field strength|29j. 

Now we go to a lattice QCD framework and perform 
numerical simulations in pure SU (2) QCD. We adopt an 
improved Iwasaki gluonic action |17|. Here we use ther- 
malized 2000 vacuum configurations at the lattice dis- 
tance a = 0.0792(2) fm. Simulation details are the same 
as in Ref.Q. 

A non-Abelian field strength F /JjV (s) is given by a 1 x 1 
plaquette variable defined by a path-ordered product of 
four non-Abelian link matrices on the lattice: 

U^{s) = exp (iF^(s)) = U° v (s) + iU^ u (s)a a . 

The unit vector in color space is 



(6) 
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FIG. 1: E electric field profiles, r is a distance perpendicular 
to the QQ axis and W(R x T = 5 x 5) is used. The solid line 
denotes the best exponential fit. 
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FIG. 2: Monopole current distributions. 
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FIG. 3: Components of monopoles around a static quark pair. 



and the Abelian-like field strength is written similarly as 
in Eq.©, i.e., /^(s) = £„n^00f 

Let us try to measure, without any gauge-fixing, elec- 
tric and magnetic flux distributions by evaluating correla- 
tions of Wilson loops and the Abelian-like field strengths 
located in the perpendicular direction to the Wilson-loop 
plane. 

We define a gauge-invariant lattice 'monopole' in the 
same way as in Eq.©: 
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i0&vfap{s + A), 



(7) 



which satisfies A^fc^(s) = 0. A M (AJJ is a lattice forward 
(backward) derivative. 

First we show in Fig^ electric field profiles around a 
quark pair. Only the z-component of the electric field 
is non- vanishing and squeezed. The profiles are studied 
mainly on a perpendicular plane at the midpoint between 
the two quarks. Note that electric fields perpendicular 
to the QQ axis are found to be negligible. The solid line 
denotes the best exponential fit. 

Now let us study the violation of the Bianchi identity 
with respect to the Abelian-like field strength Eq.©. 
The Coulombic electric field coming from the static 
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FIG. 4: The azimuthal component of the Abelian monopoles 
and the magnetic displacement current around a static quark 
pair. 
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FIG. 5: The azimuthal component of the Abelian monopoles 
and the magnetic displacement current around a static quark 
pair in the MA gauge. 



source is written in the lowest perturbation theory in 
terms of the gradient of a scalar potential. Hence it 
does not contribute to the curl of the electric field nor to 
the magnetic field in the above Bianchi identity Eq.(0}. 
The dual Meissner effect says that the squeezing of the 
electric flux occurs due to cancellation of the Coulombic 
electric fields and those from solenoidal magnetic cur- 
rents. It is very interesting to see from FigO that in 
this gauge-invariant case, the gauge-invariant 'monopole' 
Eq.Q) plays the role of the solenoidal current. This is 
qualitatively similar to the monopole behaviors in the 
MA gauge B E El • 

Let us see also the r dependence of the 'monopole' 
distribution shown in Fig [3] All r and z components 
of each term in Eq.|0J are almost vanishing consistently 
with Fig[3 The magnetic displacement current A4B are 
found to be ne glig ible numerically as similarly as in the 
MA gauge |l8L Il9i 12^. We have measured also correla- 
tions between Wilson loops and electric currents defined 
by j 11 — &-vfnu- They are found to be negligible. 




FIG. 6: The correlation between the Wilson loop and the 
squared monopole density. W(R — 5,T — 5) is used. The 
solid line denotes the best exponential fit. 
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FIG. 7: Behaviors of A r A' r E z + (l/r)A r E z . The solid line 
denotes the best exponential fit. 



For the sake of comparison, we also discuss the MA 
gauge where ^2 S Tr[U )Jj (s)a3UJ l (s)a3} is maximized. In 
the MA gauge, an Abelian link variable 0^f A (s) is defined 
by a phase of the diagonal part of a non-Abelian link 
field after the gauge-fixing. An Abelian field strength 
0ff A (s) in MA gauge is defined as 6^ A (s) = 9ff A (s) + 
0f IA (s + A) - &ff A (s +£>)- 0^ IA (s). In this case, we use 
only the third component W 3, of the non-Abelian Wilson 
loop as a source. We show in Fig0]and Fig|5] azimuthal 
components of all three terms of Eq.Q in this gauge- 
invariant case and in the MA-gauge case. It is interesting 
that the peak positions of gauge-invariant and k^ IA 
in the MA gauge look similar around 0.15 fm, although 
the height and the shapes seem different. Note that in 
the MA gauge, W 3 are used as a source of a quark and 
an anti-quark. 

Let us next try to fix the type of the vacuum of pure 
SU(2) QCD. The penetration length A is determined by 
making an exponential fit to the electric field flux for 
large r regions. The best fitting curve is also plotted in 
Fig. ^ from which we fix the penetration length A [3(J • 
Next we derive the coherence length £. In a separate pa- 
per [2^, we have shown that the coherence length can 
be fixed by a measurement of the squared monopole den- 
sity k 2 (s) = Ylu=i kfii 3 ) around the QQ pair. The same 
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situation is expected in this gauge-invariant case. The 
correlation between the Wilson loop and the squared 
'monopole' density is plotted in Fig. From the ex- 
ponential fit, we may fix the coherence length. We get 
A = 0.085(5) fm and £,/V2 = 0.10(1) fm. Although the R 
and T dependences of both lengths are not studied yet, 
the value of the coherence length looks almost the same 
as that of the penetration length within the error bars. 
Hence if the same situations will continue for larger R 
in the confining string region, the type of the vacuum is 
fixed to be near the border between the type 1 and the 
type 2. This is consistent with the result of our previ- 
ous paper |2j| and the results (see [24| and references 
therein) obtained in the MA gauge and Landau gauge. 
Also similar results were obtained in SU(3) dual QCD 
in the continuum |25| . It should be stressed, however, 
that the present result is obtained in the gauge-invariant 
framework. 

Next we study an equation similar to the dual London 
equation. From Eq.Q), we get 

4ttV x k = V(V -E)- AE + V x {d A B). 

Evaluating correlations between a Wilson loop and each 
term in the righthand side, we find only z components are 

relevant: 4tt(V x k) z A r A' r E z - (l/r)A r E z . FigO 

shows the data with the use of 4000 thcrmalizcd configu- 
rations. The exponential fit can be done for large r > 3a, 
although the errors are large. The correlation length 



0.09(2) fm is compatible with the penetration length. We 
also see 4?r(V x k) z ~ m 2 E z with mT 1 = 0.10(4) fm [24| 
which is also compatible with the penetration length. 

Finally two comments are in order: (1) The Abelian- 
like field strength Eq.Q is reduced to an Abelian one 
if off-diagonal components are negligible or non- Abelian 
lxl plaquettes are well approximated by Abelian ones. 
The former occurs in the MA gauge, whereas the latter 
does in the maximally Abelian Wilson loop gauge j2|j 
where almost the same fine results as in the MA gauge 
are observed. In these cases, we may adopt ~ 1, since 
only the diagonal components are dominant. Hence the 
present gauge-invariant results could explain why only 
restricted Abelian projection schemes like the MA gauge 
look nice among infinite possible candidates. 
(2) It is very interesting that the non- Abelian action is 
written by which can be decomposed into 'monopole' 
and 'electric-current' parts with the help of the Hodge 
decomposition. Using the Michael action sum-rule |27j . 
we can evaluate the 'monopole' contribution to the static 
potential. The work is in progress. 

The numerical simulations of this work were done us- 
ing RSCC computer clusters in RIKEN. The authors 
would like to thank RIKEN for their support of com- 
puter facilities. T.S. is supported by JSPS Grant-in-Aid 
for Scientific Research on Priority Areas 13135210 and 
(B) 15340073. 
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